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Abstract. Idealizing matter as a pressureless fluid and rep- 
resenting its motion by a peculiar-velocity fleld superimposed 
on a homogeneous and isotropic Hubble expansion, we apply 
(Lagrangian) spatial averaging on an arbitrary domain D to 
the (nonlinear) equations of Newtonian cosmology and derive 
an exact, general equation for the evolution of the (domain 
dependent) scale factor a-oit'). We consider the effect of inho- 
mogeneities on the average expansion and discuss under which 
circumstances the standard description of the average motion 
in terms of Friedmann's equation holds. We find that this ef- 
fect vanishes for spatially compact models if one averages over 
the whole space. For spatially infinite inhomogeneous mod- 
els obeying the cosmological principle of large-scale isotropy 
and homogeneity, Friedmann models may provide an approx- 
imation to the average motion on the largest scales, whereas 
for hierarchical (Charlier-type) models the general expansion 
equation shows how inhomogeneities might appreciably affect 
the expansion at all scales. An averaged vorticity evolution law 
is also given. Since we employ spatial averaging, the problem of 
justifying ensemble averaging does not arise. A generalization 
of the expansion law to general relativity is straightforward for 
the case of irrotational flows and will be discussed. The effect 
may have important consequences for a variety of problems in 
large-scale structure modeling as well as for the interpretation 
of observations. 
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1. Outline of the problem 

Traditionally, cosmological models have been based on the as- 
sumption that, on a sufficiently large scale, the Universe is 
isotropic and homogeneous. As long as inhomogeneities are 
small and described as linear perturbations off a Friedmann 
model which average to zero, this picture is consistent; by 
construction rather than derivation the averaged variables are 
given by a Friedmann model. However, in the observed part of 
the Universe the matter inhomogeneities are not small (e.g., 
the r.m.s. density fluctuations are much larger than unity be- 
low the scale of clusters of galaxies), and contemporary theories 
of structure formation follow the evolution of inhomogeneities 
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into the nonlinear regime, mostly in the framework of Newto- 
nian cosmology, either with analytical approximations or nu- 
merical N-body simulations. Again, both methods to simulate 
nonlinear structure formation are constructed in such a way 
that the average flow obeys the homogeneous and isotropic 
Friedmann models. 

Basic, nontrivial questions which are usually not even raised 
in the traditional approach are, whether and how a general 
anisotropic and inhomogeneous solution can be split into an 
isotropic and homogeneous average fleld and deviations thereof 
and whether, if so, the average variables obey Friedmann's 
equation. Because of the nonlinearity of Einstein's as well as 
Newton's laws for gravitationally interacting systems, the an- 
swers to the last questions are not obvious, as has been em- 
phasized particularly by Ellis (1984). 

In contemporary models of structure formation the time- 
and length-scales as well as the amplitude of the initial fluc- 
tuations are expressed in terms of quantities of their assumed 
Friedmann backgrounds such as a{t), H{t), Q{t), etc. . Obser- 
vational data for these parameters are interpreted accordingly. 
This procedure which rules a variety of problems in structure 
formation theories (like the question whether nonbaryonic dark 
matter is needed to explain present day structure) excludes 
by assumption all inhomogeneous models which do not obey 
Friedmann's equations if averaged on some large scale, even 
those which are kinematically isotropic and homogeneous. 

In general relativity, the problem of averaging is very in- 
volved because (i) in a generic spacetime there are no pre- 
ferred time-slices one could average over. (It should even be 
difficult to recognize a Friedmann model as such if one were 
given a complicated slice of it.) (ii) the metric is a dynamical 
variable entering the field equations nonlinearly, and it is dif- 
ficult to average (or instead deform, see Carfora and Marzuoli 
1988, Carfora et al. 1990) it, (iii) a gauge problem arises in 
relating the "true" and the averaged metric. Indications that 
the "backreaction" of inhomogeneities on the global expan- 
sion may have important consequences for the structure forma- 
tion process have been put forward by Futamase (1989, 1996), 
Bildhauer (1990) and Bildhauer & Futamase (1991a, b) who 
have studied the "backreaction" effect quantitatively based on 
perturbative calculations. In particular, they found that the 
expansion is accelerated significantly where inhomogeneities 
form. Recently, renormalization group techniques in relation 
to the averaging problem have been advanced by Carfora and 
Piotrkowska (1995) who investigated in detail the connection 
between (spatial) manifold deformations and spatial averag- 
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ing at one instant, using the constraint equations of general 
relativity. 

In this paper we address the problem of the effect of non- 
linear inhomogeneities on the average expansion in Newtonian 
cosmology where also the evolution equations can be averaged 
without uncontrolled approximations. This provides straight- 
forward insight and some answers which arc indicative also 
for the GR-casc and uncovers possible limitations of current 
cosmological models. Moreover, we will show that for scalar 
quantities (like the expansion rate) the averaging procedure 
proposed in the present work carries over to the general rela- 
tivistic case, if we restrict the equations to irrotational flows. 



2. Averages in Newtonian cosmology 

According to Newtonian physics, the motion of a sclf-gravitat- 
ing, prossurcless fluid ("dust") is governed by the Euler- 
Poisson system of equations. Thus, with respect to a non- 
rotating Eulerian coordinate system^ the fields of mass den- 
sity Q{x,t) > 0, velocity v{x,t) and gravitational acceleration 
g{x, t) are required to satisfy 



dtv = —{v ■ 'V)v + g 
dtQ = -V • {qv) , 

VX£, = 0, 

V • 5 = A - At^Gq , 



(la) 
(Ife) 
(Ic) 
(Id) 



where A denotes the cosmological constant, here included for 
the sake of generality. 

It is useful to introduce the rates of expansion 6 = V ■ v, 
shear a = (aij) and rotation w = |V x d of the fluid flow via 
the decomposition* 

Vij = <7ij + -5ij9 + u!ij ; ujij — —riijku>k , i^iij] — , (2a) 

of the velocity gradient, where 6ij and rnjk denote the (Eu- 
clidean) spatial metric and its volume form (Levi-Civita ten- 
sor), respectively. In contrast to v, the tensor fields ct, 6 and 
u! are independent of the (non-rotating) coordinate system; v 
can be reconstructed from a, 9 and w up to a spatially constant 
summand (compare Appendix A) . 

Using the Lagrangian time derivative operator ^ = (...)' =: 
dt + v-V, we may replace the system (1) by the equivalent sys- 
tem (2) consisting of eq. (2a) and the transport equations^ 
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w = —-u6 -h £ • w , 

o 

e = A- inGg - + 2(a;^ - a^) , 



(26) 
(2c) 

(2d) 



for g, u> and 6, in which a and uj denote the magnitudes of 
shear and rotation, respectively: 

er = Id := {^cnjaijf^ , w = |w| =: {^uJijUHj)^^'^ . 



^ notes are given at the end of the paper. 
* a comma denotes partial differentiation with respect to Eu- 
lerian coordinates; we adopt the summation convention. 



((1) i.s equivalent to (2) as a system for g and v. In fact, if 
(la) is taken to define g in terms of v, then (Ic) 4^ (2c) and 
(Id) ^ (2d); (16) ^ (26) is obvious.) 

We wish to take spatial averages of the equations (2b,c,d). 
For this purpose, lot f^:Xi-~^x denote the mapping which 
takes initial positions X of fluid particles at time to to their 
positions x at time t; in other words, let x = ft{X) = f{X,t) 
be the field of trajectories on which the Lagrangian description 
of fiuid motion is based. Then, the volume elements at t and 
to are related by d^x = Jd^X, where J{X,t) := is the 
Jacobian determinant of /(. Therefore, 



J = J9 , 



(3) 



and the "comoving" volume V{t) =: a^{t) of a compact por- 
tion ©(t) of the fluid changes according to 



V 



^ f d'x = f d'X j= j d 

"•^ J-D(t) ^B(to) J'D^t) 



X 9 



l-D(t) J-D(to) 

which may be written 

m\ - ^ - q'^^ 
V a-D 



(4) 



Here and in the sequel, {A}^ = ^ J-p d^x A denotes the spa- 
tial average of a (spatial) tensor field A on the domain T>{i) 
occupied by the amount of fluid considered, and ar>(t) is the 
scale factor of that domain. Suppressing temporarily the index 
V for notational simplicity, we note the useful "commutation 
rule" 



{AY - {A) = {A9) - {A){9) . 
(Proof: 



(5) 



{AY 



((y-1 / d^x A) = -^{A) + V-^ J dPX {AJ + Aj); 



using (3) and (4) gives (5). 

Applied in turn to g, uj and 9, eq. (5) combined with 
(2b,c,d), respectively, leads to 



{qY 



-{Q){e) 



(w)- = (w Vi;) - (w)(e) 



(6) 
(7) 



(0)- = A - 47rG(e) + \{{9')) - {9f + 2({a;^) - {a')) , (8) 

where we have also used eq. (2a). 

In terms of the scale factor o, eq. (6) says that a^{g) =■ M 
is the constant total mass of the fluid portion considered. Eq. 
(8) may be rewritten as 

3^ + 4^G^ - A = I {{9^) - {9Y) + 2(a;^ - cr^) . (9a) 

The averaged Raychaudhuri equation (9a) can also be writ- 
ten as a standard EYiedmann equation for the "effective mass 
density" geflt, 

47rGeeff :=477G(e)-|((^-(e))')+2(a'-a;') . (96) 

Eq. (9) shows that inhomogeneities have an accelerating effect 
on the averaged volume expansion rate {9) , if the shear term 
dominates the vorticity and contraction terms in that equation. 
Eqs. (9) can also be used to discuss anisotropics: If the domain 
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■P is chosen as a cone with the observer at the vertex, (9a) 
shows that the average expansion in two such cones will be 
different if the averages involving 0, u and a differ. 

Prom eq. (7) we conclude that 



(w) = (oj • Vv) = 



(fOa) 



The converse only holds asymptotically for positive average 

expansion: 

if (w ■ Vv) = 0, then 



, . ^, - r d* <"> (*-*0) 



(106) 



for a universe which collapses on average, (w) blows up (com- 
pare Olson & Sachs 1973)^. 



the Hubble flow is irrotational, we have w = x m; then, 
using Green's formula. 



(aj)x) = iaj,^ / dS X u , 
^ Jd-D(t) 



(15) 



((w)x))' + 2Ji'(w)D = (w ■ Vw)i> — (w)i>aj,^ / dS ■ u . 

Jd-D{t) 

(16) 

(Olson & Sachs (1973) have derived an evolution equation for 
the ensemble average of assuming homogeneous-isotropic 
turbulence.) 



The equations (6)-(9) hold for any part of any pressureless 
fluid. 

We want to apply the foregoing results to a large, typical 
part of a cosmological model. For this purpose we split the ve- 
locity field into a Hubble flow and a peculiar velocity u. Recall 
that a Hubble flow can be characterized as having vanishing 
rates of shear and rotation and a spatially constant rate of 
expansion 3H{t) (cf. eq. (2a)). Therefore, 



Vij = H{t)Sij + Ui, 
hence 

e = 3H + V-u . 
We may choose* 

or, equivalently, 
(V ■ = . 



(11a) 



(116) 



(12o) 



(126) 



H{t) and u then depend, of course, on the choice of the aver- 
aging region V. 

From eqs. (lla,b) and (2a) we get by a little computation an 
expression for the quadratic principal invariant of the velocity 
gradient {vij): 



= 6H^ +4:H'V -u + V ■ {uV ■ u-u-Vu) . 

We use it and eq. (12a) to rewrite the averaged Raychaudhuri 
equation (9a) in the form 



a-D 



a-j^ / dS ■ {uV ■ u — u ■ Vw) , 

Jd-D(t) 



(14) 

where we have reinserted the index T> to exhibit the dependence 
on the comoving domain chosen, and we have applied Gaufi's 
theorem to transform the volume integral in the average to a 
surface integral over the boundary dT> of the domain. 

We note that the averaged Helmholtz vorticity equation (7) 
can similarly be brought into a more transparent form. Since 



We emphasize that, given a solution (p, v) of the basic equa- 
tions (1), the terms in eqs. (12)-(16) depend on the choice of 
an "initial" domain ©(to) and a choice of a Hubble function 
H{t). 

So far, all equations refer to a single non-rotating, rectan- 
gular Eulerian coordinate system a". Such a description holds 
globally, if space is assumed to be the standard Euclidean M^. 
However, we wish to consider also spatially compact Newto- 
nian models as defined and analyzed in (Brauer et al. 1994). 
As shown there, space can then be considered without loss of 
generality as a torus T^, which cannot be covered by a single 
coordinate system. It is important (especially for Subsection 
3.1) that nevertheless the eqs. (14)-(16) remain valid for such 
models as intrinsic, coordinate free relations for any compact, 
oricntable domain 2?(to) and its evolution 2?(t), provided the 
models admit a Hubble fiow. The reason is that, although the 
inhomogeneous fluid motion (locally given by v) as well as the 
Hubble flow cannot be described globally by smooth vector 
fields on the torus, the peculiar -velocity u, being the relative 
velocity of those two motions, is a smooth, global 3— vector 
field on (see Appendix A for a detailed discussion of the 
kinematics on toroidal models). 



3. Interpretation for different cosmologies 

Let us now assume that the Hubble flow represents the mean 
motion on some domain T> and, correspondingly, that the 
peculiar velocity field u represents the inhomogeneous fiow 
with {uij)T> = which, by eq. (11a), is equivalent to eq. (4) 
(13) and {aij)T> = 0, {uJij)T> = 0. Eq. (14) then shows how this field 



u determines the deviation of the expansion, on the scale ajy, 
from Friedmann's law. The "perturbing terms" are quadratic 
in u and are in general nonzero (even for isotropic expansion); 
a similar remark applies to eq. (16). 

We now consider three theoretically possible cosmologies. 



3.1. Spatially compact cosmologies 

These models with a closed (i.e., compact without boundary) 
3-space admit either none or exactly one (global) Hubble flow, 
for the existence of such a flow means that the toroidal space 
expands self-similarly without rotation (see Appendix A). In 
this latter case we may choose the compact domain 'D{t) to be 
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the whole space. Then, 'D{t) has no boundary, and we obtain 
from eqs. (11a) and (14)-(16): 

3«+47rG^-A = , 
a a-^ 

{aij) = ; (u;) = , (oj • Vu) = , (17a) 
and, from (9a) and (14), 

l{{e-{9)f) + {u,') = {a'} . (176) 

(Here we have used that u is globally well-defined.) 

For inhomogeneous models based on such Priedmann back- 
grounds Brauer et al. (loc.cit.) have established some exact 
qualitative results, and in (Ehlers & Bucliert 1996) we develop 
a perturbation scheme the terms of which are well-defined to 
arbitrary order and are uniquely determined by initial data; 
that scheme is also useful for numerical work. 

Toroidal models with a background Hubble flow are, in ef- 
fect, used in all current structure formation simulations, since 
they employ periodic boundary conditions for the density and 
peculiar-velocity fields on some large scale. As we have just 
shown, such (Newtonian) models can always be regarded as fi- 
nite perturbations of toroidal Friedmann models. The pertur- 
bations, whether small or large, cannot influence the overall 
expansion. 

3.2. Cosmologies based on the principle of 'large-scale 
homogeneity and isotropy' 

If we take space to be infinite Euclidean R**, we can assume 
the boundary term in (14) to be small compared to This 
may be the case if the inhomogeneities axe substantially smaller 
than the size a-p of the averaging domain and if their peculiar- 
velocities are small on the averaging scale (|tx| << a-oH); it 
corresponds to the cosmological principle of 'statistical homo- 
geneity and isotropy'. Then, the average motion may be ap- 
proximately given by a Friedmann model on a scale which is 
considerably larger than the largest existing inhomogeneities. 
This is the generally held view. 

3.3. Hierarchical (Charlier-type) cosmologies 

Observational evidence for structures on scales of hundreds of 
Megaparsecs indicates that the "closure scale" - the circum- 
ference of the Universe if it is closed - has certainly not been 
reached by observed samples of the Universe. Power spectrum 
estimates of the density contrast show a negative slope close to 
the limit of the deepest available surveys. However, the CORE 
microwave background measurement suggests that the spec- 
trum bends over to positive slopes on larger scales, supporting 
the assumption of large-scale homogeneity. Thus, if the CORE 
measurement is indeed a detection of primordial density fiuc- 
tuations, then (up to the cosmic variance of this measurement) 
we may use Friedmann models as approximate models for the 
average fiow on the largest scales. 

The expansion law (14), however, indicates that Fried- 
marm's law may be appreciably modified by peculiar-velocity 
terms; for large-scale homogeneity and isotropy the sum of 
these terms will approach a limiting value which may or may 
not be negligible, in contrast to the spatially compact case. 



Only in the case where this limiting value is negligible, global 
values of, e.g., the Hubble constant and the density parame- 
ter are related to the Friedmann parameters on that "asymp- 
totic" scale. If such a scale doesn't exist, the average flow will 
evolve into an anisotropic and (after development of multi- 
streaming) rotational flow, in spite of isotropic and irrotational 
initial data. The extreme opposite, a (globally) hierarchical cos- 
mology, where the spectrum of fluctuations continues to rise 
on large scales, could only be treated by including the source 
terms. 

3.4. Supplementary remarks 

Eqs. (9a) and (14) imply that the expression 

Q:=l{{e-{e}f) + 2{cv'}-2{a') (18) 
is a divergence. From this we conclude: 

i) If the shear vanishes in a toroidal model which admits a back- 
ground Hubble flow, the model is homogeneous and isotropic; 

in other words: 

ii) In order to contain any inhomogeneities at all, compact 
models admitting a background Hubble flow must have nonva- 
nishing shear. This statement applies to the shear scalar, the 
shear tensor as well as to the average (a^). 

iii) If, on a// sufficiently large scales - e.g., because of the struc- 
ture formation process - the expression Q in nonzero, then the 
model either is not compact or does not admit a background 
Hubble flow. 

iv) Even if the expression Q is zero on some large scale A, 
then still 7^ and, in general, ^ {0}a, provided inho- 
mogeneities are present. 

These four statements are equivalent; they only emphasize dif- 
ferent aspects of i). To prove i), note that if <t = and the 

model is compact, the integral over the nonnegative function 
|((6l - {e)f) + 2{oj'^) varfishes; hence 9 = (6*) and {uj) = 0, 
therefore, v itself represents a Hubble flow, q.e.d. 

Note that Remarks ii) and iv) concerning the shear follow, 
because the average over tr^ is performed over non-negative 
terms and can only vanish, if the shear scalar vanishes point- 
wise and, hence, if each component of the shear teusor itself is 
zero. Turning this argument around, nonvanishing components 
of the shear tensor result in a nonzero average of on the scale 
A. Since this average does not vanish, the average fluctuation 
of the expansion does neither, unless the average vorticity ex- 
actly compensates the average shear (see also Yodzis 1974). 
In general we neither expect the summands in (18) to vanish 
individually, nor to decrease by going to larger volumes^ . This 
indicates that the "conspiracy assumption" Q = on some 
large scale (if the Universe is not genuinely compact) must be 
considered a strong restriction of generality. The existence of 
a Hubble flow (as a reference flow) does not imply that the 
average flow is a Hubble flow. 

As we have indicated in the introduction, the averaging 
problem in general relativity is very involved. However, for 
scalar quantities like the expansion rate 6 we may derive an 
expansion law in full analogy to the Newtonian case. For this 
purpose we may foliate the spacetime into a family of space-like 
hypersurfaces with spatial metric gij which are flow-orthogonal 
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(this is only possible, if the flow is irrotational). Introducing 
Gaussian (normal) coordinates we can define spatial averages 
(e.g. the average of 9) in a spatial domain V with volume 
V = dPX^, g := det(gij), as follows: 

'■=V J ^^'^'^ • (19) 

This definition agrees with equation (4), if J is replaced by 

y/g. (Indeed, if the metric is written as the quadratic form 
9ij ~ v'lv'^j involving the one forms if = r/^^dX^, then = 
det(77") =: J. ,/ is identical to the Jacobian determinant used in 
the present work, if r;" = Vo/.) Since Raychaudhuri's equation 
is the same in GR we conclude that the expansion law (9) is 
also valid in general relativity. 

Yodzis (1974) has shown (Theorem 2) that, for general- 
rolativistic, irrotational dust models for which the hypersur- 
faces orthogonal to the dust world lines are closed and ori- 
entable, equation (9) holds if the averages are performed over 
the whole space. Our statement is more general in the sense 
that equation (9) holds in general relativity for the same as- 
sumptions, but for arbitrarily chosen spatial domains. 

An important difference to the Newtonian treatment, be- 
sides spatial curvature, arises due to the fact that it may not be 
in general possible to represent the term (18) as a divergence 
in GR. We stress that this would imply a strong challenge for 
the standard cosmologies, since we can no longer argue, except 
for non-generic situations, that there exist cases in which the 
average obeys Friedmann's law. Even more, we don't expect 
the previously discussed arguments (after eq. (18)) to hold, 
since the valid theory on the large scales under consideration 
is general relativity. 
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4. Notes 

1) For spatially unbounded mass distributions there are no 
inertial coordinate systems related by Galilean transfor- 
mations. Instead, the preferred coordinate systems {t,x'^) 
are related by the more general transformations a;" = 
R\x^ -f (t), where R], denotes a constant rotation ma- 
trix and d"" {t) represents a translation depending arbitrar- 
ily on time t. The basic equations (1) are covariant with 
respect to these transformations, provided g, v and g are 
transformed in the obvious way. These coordinate systems 
are called (dynamically) non-rotating, since with respect to 
them, no Coriolis forces occur (Heckmann and Schiicking 
1955,56). 

2) We call eqs. (2b,c,d) transport equations since they describe 
how Q, LJ, 9 change along a fluid trajectory. In Lagrangian 
coordinates (t, X) eqs. (2) are ordinary differential equar 
tions with X as parameters. 

3) The assumption of homogeneous-isotropic turbulence em- 
ployed by Olson & Sachs (1973) imposes from the start 



a strong restriction which excludes the relevant terms dis- 
cussed here. Also, they use ensemble averaging instead of 
spatial averaging. Of course, statistical statements for spa- 
tial averages can be investigated by averaging over statis- 
tical ensembles of spatial domains T) centered at random 
points in space. 

4) Given Vi and H(t), (11a) defines u up to a spatially con- 
stant term. One can fix u uniquely by requiring either 
(m) = or {qu) = 0. a weaker assumption than eq. (12b) 
is to require (V • u) =0 only on the largest scales. In 
this case H{t) could be defined in terms of a Friedmann 
solution, whereas Hxi '■= could be interpreted as the 
value of Bubble's constant as measured on smaller scales 
(after averaging over statistical ensembles of spatial do- 
mains). The inhomogeneous term on the r.h.s. of eq. (14) 
would then contain additionally a surface integral of the 
flux of u through the boundary &D. Since a-p is deflned via 
the volume, it could be interpreted also as an "effective" 
(D-dependent) scale factor of a possibly anisotropically ex- 
panding domain (compare also Appendix B and, for further 
discussion, Buchert (1996)). 

5) Imagine we divide a large portion of the Universe "D with 
volume V{t) into A'^ subdomains Vi with volumes Vi{t). 
Then, e.g., (c^)x> = If the subdomains are 
'typical', then we conclude (cr^)xi = (it^)d^, i.e., the aver- 
age value of (T^ attained in the subdomains is "frozen" and 
cannot become smaller by averaging over larger domains. 
This statement also applies to any positive semi-definite 
quantity such as lj^ and {9 — (9))^. 
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Appendix A: Kinematics on toroidal spacetimes 

Consider any "Newtonian" toroidal spacetime § = x R . 
The "global motion" of the space T"* of § is given by a time- 
dependent triad e; {t) of spatial vectors generating those trans- 
lations n''ei{t) which map each point of the torus onto itself 
at time t, being arbitrary integers. The rates of change e^, 
measured with respect to a non-rotating orthonormal basis, 
uniquely determine a matrix Hij{t) via ej = HijBj. 

With respect to any particle as origin, any point of has 
infinitely many position vectors x + n^ei{t) at time t. Hence, 
any kinematically possible motion covering the universe § can 
be represented by a velocity field v{x^ t) on the covering space- 
time X R of S which obeys the torus condition 

v{x + n^ei{t),t)=v{x,t)+n^ei{t) . (T) 

It relates the velocities associated with the position vectors of 
any one particle. The diff'crcnce of two such velocity fields, but 
not these fields themselves, define global spatial vector fields on 
the torus, as (T) shows. This is illustrated in the figure below: 




If V obeys (T), then Vij{x + n^ei{t),t) = Vij{x,t), hence 
(7ij, ujij and are tensor fields on T^. 

A velocity field V obeys the cosmological homogeneity prin- 
ciple iff, for all y and x, 

Viy,t)-V{x,t) = A(t){y-x) (C) 

(Heckmann & Schiicking 1955, 1956). This property is compat- 
ible with (T) iS A = H, as follows by putting y = x + n^ei{t) in 
equation (C) and using (T). This shows: Except for a spatially 
constant additive contribution, there exists exactly one homo- 
geneous motion on a toroidal spacetime; a representative of its 
velocity field is determined by the matrix Hij defined above, 
as V{x, t) = H{t)x. Its shear Sj^, rotation Q,ij and expansion 
G are spatially constant and given by the decomposition 

o 

If V is any inhomogeneous fiow on S and V + A(t) the general 
homogeneous one, then a unique peculiar-velocity field u is 
defined hy u := v — V , with {u)j3 — 0. Since w is a global 
vector field on T'"*, {uij)j3 = 0, hence 

{u}ij)T3 = Oij , {aij)r3 = Eij , (9)^3 = 6 . 

The homogeneous motion of a toroidal model is a Hubble flow if 
and only if Eij = and ilij = or, equivalently, Hij = ^QSij; 

the last condition is in turn equivalent to ei{t) = ^^ei(*o), 

if we put e = Sf. 



Thus, a Hubble flow exists and is then uniquely determined on 
a toroidal model if and only if the torus expands self-similarly 
without rotation. 

Appendix B: General expansion law with global shear 
and vorticity 

We start from Raychaudhuri's equation, averaged on a simply 
connected domain (eq. (9a)), 

3^ + 47rG^ - A = ^ {{9')^ - (9)1) + 2{J' - a^)^, . 

Qi'D 

By introducing on V arbitrary fields of expansion G = G(t), 

shear = Eij(t) and vorticity Q.ij = 0.ij(t) (which are not 
necessarily average values, but merely define a time-dependent 
standard of reference, e.g., a homogeneous solution of the ba- 
sic system of equations), we may define a linear "background 
velocity field" V by Vi := HijXj with 

Vij = Eij -|- —QSij -\- ilij =: Hij . (B-l) 

Also, the inhomogeneous deviations from these reference values 
are introduced, 

9 = Q + 9 ; Oij = Yjij + aij ; tOij — Q,ij -\- uJij . {B.2a, b, c) 
Equation (9a) may then be cast into the form 

3^+47rG^-A = 
a-D aj, 

2 (Sl'^ — E'^) + 2 {^iij{uJij)T> — Ejj(6'ij)Ei) 
+ li{e'}T,-{9}l)+2{u,^-a^}T,. 

(5.3) 

Using (13) to express deviations from homogeneity in terms 
of the peculiar-velocity gradient (uij) with H^ = Hij Hij = 
2(E^ -|- n^) -|- |G^, we arrive at the (most general) expansion 
law which holds for general inhomogeneities relative to a non- 
isotropic and rotational "background velocity field" : 

a-D aj, 
2 (q'^ - E^) + 2 (fl^j{Lb^j)T} - E,,((T,j)i,) 

-KV-[«(V-w)-(wV)«])D-|{V-M)i . (B.4) 

Instead of considering inhomogeneities relative to a (D-inde- 
pendent) arbitrary reference velocity field V, we may alter- 
natively define the inhomogeneous fields in (B.2a,b,c) as de- 
viations from the average flow within the domain X> (which 
we identify with the global background velocity field). This 
corresponds to the point of view of an observer who maps a 
finite region of space and assumes that the average values of 
that region are 'typical' for the Universe. In this case all global 
variables @{t), Eij(t) and ^iij{t) are averages and depend on 
content, shape and position of the spatial domain "D: 

= {9)t> ; T,ij = {(7ij)T> ; fiij = (wij)© . (B.5o, 6, c) 
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The averages {6)t> = (V • u)t>, {d'ij)'D and {uiij)^ then vanish 
by definition, and the expansion law (B.4) simpUfies to 

3 VAttG^ — A = 

a-D a'±i 

2(n^ - E^) + (V • [m(V • It) - (m • V)u])t> ■ (B.6) 

The last term in (B.6) is, via Gaufi's theorem, a surface integral 
over the boundary of V. In case of a toroidal model we may 
choose T) to be the whole torus. Then, the background fiow is 



unique (see Appendix A) and, since there is no boundary, we 
obtain the global expansion law (in agreement with the result 
of Brauer et al. loc. cit., eq.(2.15)): 

3«+47rG4-A = 2(n2-E2) . (B.6') 
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